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We show that a bilayer graphene flake deposited above a ferromagnetic insulator can behave as a spin-filtering 
device. The ferromagnetic material induces exchange splitting in the graphene flake, and due to the Fano 
antiresonances occurring in the transmission of the graphene flake as a function of flake length and energy, it 
is possible to obtain a net spin current. This happens when an antiresonance for one spin channel coincides 
with a maximum transmission for the opposite spin. We propose these structures as a means to obtain 
spin-polarized currents and spin filters in graphene-based systems. 
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I. INTRODUCTION 

In the last years, there has been much interest in ex- 
ploring the unique properties of nanostructures for spin- 
tronic devices, which utilize the spin degree of freedom of 
the electron as the basis of their operation^. To this end, 
novel ways of generating and detecting spin-polarized 
currents have been explored. For instance, Song et al£ 
described how a spin filter might be achieved in open sys- 
tems by exploiting the Fano resonances occurring in their 
transmission characteristics. The Fano effect^ arises 
when quantum interference takes place between two com- 
peting pathways, one connecting with a continuum of en- 
ergy states and the other with discrete states. In open 
systems in which the spin degeneracy of carriers has been 
lifted, this effect might be used as an effective means to 
generate spin polarization of transmitted carriers. The 
idea is to tune the system so that a transmission reso- 
nance for one spin channel coincides with an antireso- 
nance for the opposite spin. In this way, a spin polarized 
current arises. 

Previous works have shown that graphene bilayer 
flakes exhibit Fano antiresonances in the transmission^. 
We propose to exploit these antiresonances to produce 
spin polarized currents in a graphene-based system by 
putting the graphene flake in contact with a magnetic 
insulator, such as EuO?>®. Exchange splitting is in- 
duced in the graphene flake due to the magnetic prox- 
imity effect— . This produces an opposite energy shift 
in the spin-up and down antiresonances in the con- 
ductance, yielding a spin polarization of the current. 
The feasibility of spin-polarized currents in graphene de- 
vices is important for the development of all-graphene 
electronics, which is one of the goals in the research 
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FIG. 1. Schematic view of a graphene flake deposited onto 
a graphene nanoribbon. This system can be considered as 
a bilayer flake with monolayer ribbon contacts. The bilayer 
portion of the system is placed over a ferromagnetic material 
(blue). 



focused on graphene applications 12 . Recently, it has 
been experimentally demonstrated the growth of Eu 
on graphene, thus bringing closer the possibility of the 
proposed scheme to induce spin polarized currents in 
graphene-based systems M- 

In this work we study the transport properties of a 
ferromagnetic bilayer graphene flake, showing that it can 
behave as a spin filter. We consider that the flake is con- 
tacted by monolayer nanoribbons, and the contacts are 
connected to the same layer of the flake, as shown in fig- 
ure [T] Such configuration can be achieved by placing a 
monolayer flake onto a nanoribbon. It is also possible 
to form a bilayer flake by overlapping two semi-infinite 
nanoribbons£*«; however, we concentrate here in the 
bottom-bottom geometry, given that the changes due to 
the configuration of the contacts merely produce a varia- 
tion in the position of the transmission antiresonances^. 
We employ a one-orbital tight-binding model which we 
solve analytically within the single-mode approximation. 
This solution agrees perfectly with the numerical result 
obtained by a Green function method in the one-mode 
energy range. 

Our main results are the following: 
(i) We obtain an analytical solution for the spin- 
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FIG. 2. Transmission versus energy for spin-up carriers of 
the armchair bilayer flake of width W = 5 and length L = 30 
(black thin line) and the coupled-chain model of equivalent 
length, n — 120 (blue thick line). The latter is shown only for 
the one-mode energy interval of the armchair ribbon system. 



dependent transmission through bilayer graphcne flakes. 
The comparison of this analytical result to the numer- 
ically computed transmission, obtained by a recursive 
Green function method, is excellent in the one-mode en- 
ergy range. 

(ii) The analytical expressions for the transmission allows 
us to explore thoroughly the parameter space, locating 
the most advantageous system sizes to obtain a net spin 
current. 

(iii) We have found that the maximum spin polarization 
is obtained when sharp antiresonances are produced in a 
plateau with maximum transmission. These correspond 
to quasi-localized states in the bilayer graphene flake. 
Thus, the tuning of the flake length is important to ob- 
tain a net spin current. 



II. MODEL 

As we focus on the electronic properties close to the 
Fermi energy, we employ a 7r-orbital tight-binding model, 
which gives an excellent description of the low-energy 
properties of graphene systems. We consider the in-plane 
nearest-neighbor interaction given by a single hopping 
parameter 70, which is approximately 3 eV, concentrat- 
ing on metallic armchair nanoribbons, which can play 
the role of contacts in this system. We would like to 
note that, in principle, all armchair nanoribbons show 
a gap due to edge reconstruction, as evidenced by first- 
principles calculations^ 5 -. A change of the hopping pa- 
rameter at the edges can reproduce this effect within 
the one-orbital approximatio n 15 ! 16 ; however, we choose 
all hoppings to be equal in order to solve analytically 
the problem and then compare to the numerical calcu- 



lations. Just above the gap, which diminishes with the 
ribbon width, the behaviour is as reported here. 

With respect to the bilayer flake, we focus in one type 
of stacking, the AA, taking into account only one hop- 
ping parameter "fi between atoms placed directly on top 
of each other. We assume that 71 = O.I70. Note that, 
although the AB stacking is more stable for graphite, 
the direct or AA stacking has been experimentally found 
in few-layer graphene^. Moreover, the antiresonances 
which give rise to the energy windows where spin po- 
larization takes place in these systems are present both 
in AA- or in AB-stacked bilayer graphene flakes, so the 
physics of the problem can be studied in any of the 
two geometries. Our choice of the AA stacking is mo- 
tivated by that fact that, as we show in this paper, the 
AA-stacked system can be easily mapped into a one- 
dimensional chain with two hoppings, yielding analytical 
expressions of the transmission that perfectly fit the con- 
ductance obtained numerically in the single mode regime. 
Within this description, we measure the length of the 
flake in translational unit cells, which corresponds to the 
number of 4-atom units along the flake. For the width 
we use the standard notation, giving the number of dimcr 
chains across the flake. 

The Hamiltonian of the complete system can be writ- 
ten as a sum of the left and right lead parts plus a central 
part, which constitutes the bilayer flake: 

n = H L + -H R + H F (1) 

with 

Ul,r = -70 Y C L C J>' (2) 

where c\ a (ci >a ) stands for the creation (annihilation) 
operator of an electron with spin a on site i, and the 
sum is over nearest neighbors. The Hamiltonian of 
the central conductor, i.e., the bilayer graphene flake, is 
given by 

U F = -7o Y c L,a c J>,« - 7i XX c i,°-,tt Ci >°-' d + h ' c ') 

+ Y 17 A <* C U,a £ W ( 3 ) 

In this latter expression a runs over u, d, i.e., the up- 
per and lower flake, respectively, the spin index a takes 
the values ±1, and A cx is the exchange splitting induced 
by the ferromagnetic insulator. For graphene systems on 
EuO we take the value A cx = 0.00l7cr>£. For the sake of 
simplicity in the analytical approach, we assume that the 
same ferromagnetic term is included in both layers of the 
flake. However, we have numerically computed both pos- 
sibilities (i.e., spin interaction in one or in two layers) and 
found the main difference is that the energy shift of the 
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antiresonances is twice when the ferromagnetic coupling 
affects the two layers. 

The transport properties of this structure can be ob- 
tained numerically by computing the Green function 
of the system, which yields the conductance in the 
Landauer-Kubo approximation 18 . However, as this sys- 



tem is uniform with a high symmetry, one can take ad- 
vantage of this and employ an analytical approach to 
solve the problem in the single-mode approximation. 

The tight-binding equations of motion for the ampli- 
tude of probability to find an electron in site of atom 
A(B) in a bilayer ribbon flake are given by 
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Assuming that the solutions in the y-direction are a com- 
bination of plane waves of the form e ±mq , the expressions 



(H]) are reduced to the equations of motion of two coupled 
dimer chains (see Appendix), which eventually lead to a 
spin-dependent transmission given by 
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In the above U„(x) are the Chebyshev polynomials of the 
second kind and E a ,± = E — a A ox =p 7i- 

In order to quantify the spin-dependent transport, it 
is necessary to evaluate the degree of spin polarization of 
the electric current. There are several ways to do this, 
but following Ref.— we introduce the weighted spin po- 
larization as: 



p \ T t ~ T i\ rp 



(7) 



Notice that this definition takes into account not only 
the relative fraction of one of the spins, but also the con- 
tribution of those spins to the electric current. In other 
words, we will require that not only the first term of the 
right-hand side of ([7]) to be of order of unity, but also the 
transmission probability T a (E) should be significant. 



III. RESULTS 

We begin by comparing the analytical results for 
the one-dimensional chain to the conductance of the 
graphene nanoribbon system obtained numerically em- 
ploying a decimation technique to obtain the Green func- 
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FIG. 3. Contour plot of the transmission versus energy and 
coupled-chain length n without ferromagnetic gate (A cx = 0). 



tion of the structure, as done previously for similar sys- 
tems without spin polarization^^. Given that the shift 
produced by the effective exchange is small, we just 
present the results for the spin-up channel. As the arm- 
chair ribbon length is measured in units comprised of 
four atoms, the equivalence between the one-dimensional 
chain and the armchair ribbon length is n — AN. Fig. [5] 
shows the transmission for the n = 120 one-dimensional 
chain model along with the transmission for the arm- 
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FIG. 4. (Color online) Spin-dependent transmission for the 
graphene flakes of length n = 40, 60 and 120 with A cx = 
O.OOI70. The solid black line is for spin up, and the solid red 
line represents the spin down. 



chair nanoribbon system of width W — 5 (two complete 
hexagons) and bilayer flake length N = 30. The energy 
interval is chosen to display the results for the armchair 
system with two modes contributing to the total trans- 
mission. As expected, the results for the one-dimensional 
chain transmission and the one obtained for the armchair 
nanoribbon system perfectly coincide in the whole one- 
mode energy window. Thus, in the following, we focus 
on the analytic results in the single-mode approximation. 

The analytical solution of the one-dimensional coupled 
chain allows us for exploring thoroughly the dependence 
of the transmission coefficient on the flake length n as 
a function of energy. As the exchange energy is very 
small, producing a shift of a few meV in the antireso- 
nances or any other features of the conductance, it is 
sufficient to study the transmission for zero exchange en- 
ergy in order to identify the situations more convenient 
for spin filtering. Figure [3] presents a contour plot of the 
transmission for the one-dimensional coupled chains as a 
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FIG. 5. Spin polarization versus energy for different flake 
lengths, with A cx = O.OOI70. 



function of energy and flake size for A cx = 0. Within this 
single-mode approximation, our results present an excel- 
lent agreement to those obtained within the continuum 
Dirac model for graphene^. 

From this figure[3]it is clear that for flake length n — 30 
there is a minimum in the transmission, and that the 
sharpest antiresonances occur for n = 60, slightly above 
120 and so on, where spin polarization will be maximized. 
The n = 40 has a noticeable oscillation of the transmis- 
sion, but does not reach the maximum value of 1. To see 
more clearly this effect, we choose three lengths for the 
flake, namely, 40, 60, and 120, and plot the transmission 
versus energy for the system placed over a ferromagnet, 
as described above. Figure 2] shows the transmission for 
spin-up and down electrons in these three cases. We see 
that the n = 40 case has a smoother oscillation that does 
not rise to transmission 1; however, both n — 60 and 
n = 120 reach the maximum value of the conductance 
and have very sharp antiresonances, as anticipated from 
the contour plot. 



5 



down 




FIG. 6. (Color online) Contour plot of weighted spin down po- 
larization for the bilayer graphene flake with A cx = O.OOI70, 
as a function of the Fermi energy and flake length. 



For these three flake lengths, we compute the spin po- 
larization as indicated in equation (J7J- Even for n — 40, 
where the transmission oscillates smoothly and the min- 
imum values for opposite spins are barely distinguish- 
able, there is a non-negligible value of the polarization, 
depicted in the top panel of figure [5j However, there is 
not appreciable spin filtering, given that the two peaks 
overlap in energy. For the other flake lengths, namely, 
n = 60 and n = 120, the antiresonances in the transmis- 
sion are very sharp and well-separated in energy. This 
yields a maximum polarization, shown in the center and 
bottom panels of figure [SJ 

The analytical solution allows us to explore thoroughly 
the system sizes and energy ranges, with the aim to iden- 
tify the optimal cases for spin filtering. Figure [6] is a 
contour plot of the weighted spin down polarization as 
a function of energy and bilayer flake length. As the 
spin-up and down antiresonances are very close in en- 
ergy, this plot is sufficient to identify the cases with max- 
imal weighted polarization. The maxima of the polariza- 
tion contour plot follow lines that correspond to quasi- 
localized states that give rise to sharp antiresonances. 

It is enlightening to compare this contour plot to the 
one depicting the transmission (figure [3]), where we can 
observe vertical wide bands of maximum transmission 
(in red) that shrink for some values of energy. The lines 
drawn by the polarization maxima go precisely through 
the points where the maximal transmission regions nar- 
row dramatically into a marked minimum. 

In the literature, it is has been showed that the Fano 
effect is robust under the several conditions, such as, the 
presence of impurities, electron-electron interaction, ex- 
ternal fields, etc X. In this sense, we believe that our 
model captures the essence of this phenomenon, which 



allow us to propose this spin filter device. 

IV. SUMMARY 

We have studied the transport properties of a bilayer 
graphene flake placed in close proximity to a ferromag- 
netic insulator. This produces an effective Zeeman split- 
ting that can be exploited to obtain a spin filtering effect. 
We have shown that the tight-binding model for the fer- 
romagnetic flake can be mapped into a one-dimensional 
chain that can be analytically solved in the single mode 
approximation. The analytical solution allows for a thor- 
ough exploration of the system parameters and the sub- 
sequent identification of the most convenient system sizes 
for a maximum spin polarization. 

Appendix A: Derivation of the transmission for the coupled 
linear chains 

In this Appendix we give in more detail the procedure 
to derive the spin-polarized transmission (J5j analytically. 
We use the plane- wave form of the solutions 

V^ = ^e ±iro *, (Al) 

which, substituted in the equations of motion yield 

(E ~ tjJAft = MB];: + 2 cosqB^J + 71^ 

(E - = 7o(^ + 2 cos qAf^) + yiB% (A2) 

We restrict our calculation to the case of metallic arm- 
chair nanoribbons and lower energies, close to the Dirac 
point (q = 7r/3). For this case, (IA2|) can be mapped 
into the equations for two linear chains by replacing 

Al' a — fa R9- Q _ fa 
^-j,a — Jl,a> j 1 (7 ~ Jl+l,cn 

(E - aA ex )f^ = 7o(/r+i >CT + f?-x,*) + lift 
(E - aA cx )f^ = 7o(/&i, CT + fti,a) + lift (A3) 
Employing the symmetric and antisymmetric combi- 
nation of the solutions, = f u ± f d , the equations of 
motion become 

(E - aA cx - 7l )/+ = 7o (/+ 1)CT + f+ ha ) 

(E - aA cx + 70/^ = 7 o(W + fU,*)- ( A4 ) 
In order to obtain the transmission through the sys- 
tem we consider open boundary conditions. Plane waves 
incoming from the left with amplitude unity are reflected 
from the scattering region (the bilayer flake) with ampli- 
tude r a , and are transmitted to the right contact with 
amplitude t a : 

ft = e lkl + r<7 e- lkl , I < 

ft=t a e m ,l>N. (A5) 

This leads to a system of equations on the reflection 
and transmission amplitudes. The transmission is given 
by T a (E) = \t a \ 2 . 
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